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A Study of the Stress Wave Factor Technique 
for the Cnaracterlzatlon of Composite Materials 


(ABSTRACT) 

This study has Investigated the potential of the Stress Wave Factor 
as an NOT technique for chin composite laminates. Tite conventional SWf 
and an alternate method for quantifying the SWF were Investigated. 

Agreement between the Initial SWF number, ultrasonic C-scan, In- 
plane displacements as obtained by full field moire Interferometry, and 
the failure location have been observed. The SWF number was observed to 
the highest when measured along the fiber direction and the lowest 
when measured across the fibers. The alternate method for quantifying 
the SWF used 1 square root of the zeroth moment (/ M Q ) of the frequency 
spectrum of the received signal as a quantitative parameter. When 
[0 f 90 2 ] s , [0,t45] s and [0,90,±45] s graphite epoxy laminates were 
cyclically loaded, 1 decreased with the number of cycles. Decreases 
1i ».he amplitudes of the higher harmonics of the frequency spectrum 
occurred as a result of transverse matrix cracking In the 90° plies of 
the laminates. The decreases in / M Q were In agreement with the 
reduction In laminate stiffness due to cyclic loading. The decreases 
in / M q were larger than the reduction in the laminate stiffnesses due 
to cyclic loading. 

From this study it therefore appears that the stress wave factor 
has an excellent potential to monitor damage development In thin 
composite laminates. 
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1. INTRODUCTION 

1.1 Scope 

The use of fiber-reinforced composite materials has become 
Increasingly popular for a wide range of applications, especially for 
high performance structures. Two of the primary advantages of composite 
materials are that they have high strength to weight and high stiffness 
to weight ratios. Also, the direction dependent properties of 
composites can be utilized in a number of applications. In order to use 
these materials effectively in engineering design one needs to 
understand: 

• how these materials behave under load. 

• the effects of initial imperfections on the final properties, 

• the nature of damage development under load, 

• the sequence of damage events that lead to the final fracture, 

® and the effect of damage on the remaining strength, stiffness 

and/or life of the material. 

Nondestructive testing (NOT) and evaluation of composite materials 
have received Increasing attention in recent years to evaluate and 
predict the mechanical properties the material in both the ar- 
received condition and during service loading. Composite materials 
develop a complex state of damage when subjected to load [1,2]. Damage 
in these materials usually includes a combination of the following* 
matrix cracks, micro delamination, macro delamination, fiber splits, 
debonding, fiber breaks, etc. Because of this complex damage state 
conventional NOT methods are unable to provide a useful parameter to 
characterize this damage. There is, therefore, a strong need to develop 
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an NDT technique which Is capable of yielding a parameter which is 
indicative of the total damage in the material. The ultrasonic Stress 
Wave Factor (SWF) is one such technique. The SWF technique was 
developed by Vary and his co-workers [3], 

The SWF technique has the potential to inspect and interrogate the 
state of damage in a material. As will be seen in section 1.2, 
ultrasound beams interact with th; entire volume of the material through 
which they travel. Therefore their propagation will be affected by the 
total damage state that lies in their path. A brief review of the work 
done in recent years in the areas of ultrasonic and acousto-ultrasonic 
NDT techniques is presented in the next section. 

1.2 Background 

Ultrasonic and more recently acousto-ultrasonic NDT techniques have 
become increasingly popular to detect damage and strength variations in 
fiber-reinforced composite materials. It has been reported [4] that 
from the measurement of extensional and shear moduli elastic and 
strength properties can be empirically evaluated. It has also been 
shown that ultrasonic attenuation in composites Is sensitive to features 
(microvoids, etc.) that influence strength properties [5], 

It has been shown that stress wave propagation is pivotal during 
failure of composites [6, 7, 8, 9]. The SWF is an aroustc 'Ultrasonic NDT 
technique. A schematic of the SWF set-up Is shown in Fig. 1 [10], A 
repeated series of ultrasonic pulses are injected into the specimen by 
means of a sending transducer. Each of these pulses simulates acoustic 
emission events in the material [19,20]. A separate receiving 
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transducer on the same side of the specimen intercepts some of the 
stress wave energy that radiates from the point of injection, the 
signals arriving at the receiving transducer resemble burst type 
acoustic emission event:* [21], The degree of attenuation of the 
ultrasonic wave as it travels through the specimen is converted to a 
numerical value called the Stress Wave Factor (SWF) after appropriate 
signal conditioning. Vary et al. [3] define the SWF as the product of 
the number of times the voltage level of a single signal exceeds a set 
threshold level, the pulse repetition rate, and a predetermined length 
of time before the counter is reset. This is schematically shown in 
Fig. 2. 

SWF = C T I (1) 

where 

C = total number of counts of a single signal that exceeds a set 
threshold level. 

T = time before the counter is reset. 

~ = pulse repetition rate. 

Vary et al . [3] describe the SWF as a measure of stress wave energy 
transmission. Their hypothesis was that decreased stress wave energy 
flow corresponds to decreased fracture resistance. The SWF provides a 
means of rating the efficiency of dynamic strain energy transfer in a 
given material. If the material exhibits an efficient stress wave 
energy transfer, then it will have higher strength. That is, better 
stress wave transmission means better transmission of dynamic stress and 
load distribution. Conversely, low values of SWF would indicate places 
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where the dynamic strain energy Is likely to concentrate and promote 
fracture. Higher attenuation gives rise to a lower SWF. In other 
words, the higher the SWF, the less the attenuation of stress waves. 
Therefore, the SWF can be used as an NDT technique. A composite with 
nigh values of SWF would therefore exhibit higher strength locally 
because the stress wave dissipation of energy is enhanced by the same 
features that increase SWF. Conversely, low values of SWF would 
Indicate that strain energy is likely to concentrate near crack 
nucleation sites and Induce crack growth [10], 

It has been shown that the SWF decreases in proportion to 
fractional powers of the ultimate strength and that it may be a useful 
aid in predicting failure locations in thin graphite epoxy laminates 
[3], The SWF was found to be a sensitive indicator of composite 
strength variations that c'.ccompany various fiber orientations relative 
to the load axis. SWF was also found to be sensitive to strength 
variations in composite laminates because of variations in the cure 
pressures [11]. The SWF was also found to correlate strongly with 
interlaminar shear strength of graphite polyimide composite laminates 
[12]. Strong correlation between the SWF and an index of fracture 
toughness was reported by Hull et al. [13], Williams et al . [14] 
reported that correlations exist between drop weight impacts, residual 
strength, through thickness attenuation and the SWF for graphite epoxy 
laminates that were subjected to repeated controlled low velocity drop 
weight impacts in the Unlnate direction. Henneke et al. [15j have 
shown that SWF can predict failure location in thin E-glass epoxy 
laminates. Govada et al. [16] reported that ultrasonic attenuation and 
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th~ SWF were sensitive to damage development In metal matrix 
composites. Duke et al . [17] reported that the SWF correlated well with 
the reduction in stiffness during fatigue loading of graphite epoxy 
composite laminates. Hayford et al. [lb] have shown that delamination 
and matrix cracking in graphite epoxy specimens can be detected by 
measuring the amplitude attenuation of an ultrasonic pulse propagated 
through the laminate thickness using a pulse-echo method. 

It Is therefore clear that ultrasonic and acousto-ultrasonic 
techniques can be used to interrogate and assess damage, and predict 
mechanical properties of composite materials. In this study the 
potential for the ultrasonic Stress Wave Factor as an NOT technique to 
assess damage and predict mechanical properties of composite materials 
is investigated. The objectives of this study are listed in the next 
section. 


1.3 


Objecti ves 
The objecti ves 
Phase I: a) 

b) 

Phase II: a) 


of this study are divided into three phases. 
Determine the feasibility of the SWF technique as a 
potential NUT method for fiber-reinforced composite 
materi als. 

Check the validity of the SWF results using other 
NDT techniques. 

Investigate an alternate method (shape parameters 
of the Fourier Transform of the acousto-ultrasonic 
signal) to quantify the SWF technique and thereby 
eliminate the need for a threshold level and 
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several other Instrumentation and experimental 
parameters. 

b) Apply the alternate method to fiber-reinforced 

composite laminates subjected to cyclic loading and 
correlate the degradation of stiffness, due to 
damage, to the root mean square value of the 
frequency spectrum. 

Phase III: Develop a computer code which can calculate various 

moments of the Fourier Transform of the received 
si gnal . 
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2. EXPERIMENTS AND RESULTS 

2.1 Phase I 

In the Initial phase of the feasibility study of the SWF as an NDT 
technique, a commercial instrument, AET model 206AU*, was used, 

Fig. 3. This Instrument contains a pulser, a receiver, a transmitter- 
receiver transducer assembly, an oscilloscope to display the received 
signal, and a counter. The specimens used were E-glass epoxy composite 
laminates with various stacking sequences. 

2.1.1 Reproducibility of the SWF Number 

For making an SWF measurement, two transducers, a transmitter and a 
receiver, must be applied to the specimen. In order to mount the 
transducers to the specimen one needs to consider the pressure applied 
to the transducers, the type and amount of couplant used between the 
transducer and the specimen, and the surface finish of the specimen. 

In order to obtain a constant reading of the SWF, it was found to 
be necessary to apply at least 2 lbs. load on each transducer. The 
value of the SWF decreased slightly with Increase in load up to about 
20 lbs. on each transducer wherefore it remained constant. For loads 
greater than 20 lbs. on each transducer, some surface damage was noticed 
on the specimen. A value of lb lbs, on each transducer was chosen, as 


Acoustic Emission Technology Corporation 
1812J Tribute Road 
Sacramento, CA 9581b 
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this *alue did not damage the specimen and the SWF number was 
sufficiency high, remained constant, and was easily reproduced. 

The various couplants used, in order of increasing viscosity, were 
Nonaq stopcock grease, U1tragel couplant, Dow-Cornlng silicone grease, 
and Panametrics shear wave couplant. Panametric couplant, which } as the 
highest viscosity, provided the highest SWF for glass epoxy specimens. 
For all the subsequent work in this study, Panametrics couplant was 
used. The amount of couplant did not significantly affect the SWF 
values or its reproducibi lity. 

The surface irregularities on the specimen due to the scrim cloth 
used in the processing procedure did not significantly affect the SWF 
values when compared with those of a polished specimen. 

Some general observations regarding the instrumentation must be 
noted: for lower repetition rates of the Input pulse there is a smaller 

scatter in the SWF values. Also, smaller scatter Is obtained at higher 
threshold levels and narrower gate widths. Figure 4 is the schematic 
representation of the interrelationship between signal, threshold, and 
gate w^dth. If the threshold is set too high and/or the gate too 
narrow, some of the useful excursions in the signal may not be Included 
in the SWF number. On the other hand, if the threshold is too low, the 
counter will detect some of the noise associated with the signal and if 
the gate is too wide it will cause the counter to detect some higher 
level noise at the tail end of the signal. Therefore a trade-off must 
be made. To obtain meaningful SWF values one must observe the shape of 
the signal and then choose the instrumentation parameters accordingly. 
The parameters chosen must yield reproducible SWF values. After 
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systematic testing and evaluation, the following Instrumentation and 
experimental parameters were found to provide reproducible SWF values to 
within +10 % for glass epoxy laminates. 


Threshold level: 1.00 volt 

Threshold mode: auto (eliminates some of the noise from the 

counting circuits 


Gate: 

Trigger rate: 
Gain: 

Sweep rate: 

Couplant: 

Pressure on each 
transducer: 


4 units 
0.5 kHz 
4b dB 

62.5 usec/dl vision 
Panametrics shear couplant 

15 lbs. 


Distance between the 
transducers: 1.5 inches 


Several instrumentation difficulties had to be overcome before the 
SWF measurements were made to be reproducible [15], The difficulties 
will not be further reported here. 

2.1.2 Correlation of SWF with Failure Location 

A large number of tensile specimens were tested for the Initial SWF 
number. These specimens were then loaded in quasi-static tension to 
failure. In general, the outcome of this testing showed that there is a 
correlation between the initial SWF number and the final failure 
location in the composite laminates. To determine this correlation, one 
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must exercise extreme care when selecting the Instrumentation and 
experimental parameters to measure the SWF*. The parameters listed in 
section 2.1.1 were used for this test program. 

The specimens tested were E-glass epoxy composite laminates. The 
dimensions of the specimens were 1.0" by 8.0" long. Several different 
stacking sequences were investigated. Measurements of the SWF were made 
along the length of the specimen at 0.25" intervals. A schematic of the 
measurement pattern Is shown in Fig. 6. Panametrics shear wave couplant 
was used. A load of 30 lbs. was applied on the transducer assembly by 
means of a spring loading frame, Fig. 3. At each location on the 
specimen, ten individual measurements of the SWF were made. After each 
measurement the transducers were removed from the specimen, the couplant 
was wiped off, fresh couplant was applied, and a new measurement at the 
same location was made. After the measurements were made, an average 
value of the SWF was determined for each location on the specimen. The 
variation of these values was within 10X. The SWF values for the 
various positions on the specimen were plotted. Figure 7 shows a 
typical normalized SWF versus position plot for a E-glasr 

laminate. Normalizing was done with respect to the maximum SWF value. 
Because of the nature of the measurement of the SWF, there is an overlap 
of values. To eliminate any shift in the data because of this overlap, 
the SWF values can be averaged. A comparison of an averaged normalized 
SWF plot and an unaveraged normalized SWF plot is shown in Fig. 8. The 
distance between the transducers is 1.5". Each measurement interval is 
0.25", The averaging is done by taking the sum of seven numbers, three 
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on either side of the location one wants to calculate the average SWF 
value, and take the mean. 

The specimens were then loaded In quasi-static tension until they 
fractured. The loading was done on an MTS servo-hydraul 1c , closed-loop 
testing frame. A schematic of each of the specimens Indicating the 
fracture, by a shaded area, was then put beneath the appropriate 
normalized SWF versus position plot. This would allow one to check for 
correlation between the lowest value of the SWF and the fracture 
location on the specimen. Figure 9 shows one such plot for a [0 * 91 ) 3 3 S 
E-glass laminate. From Fig. 9 one can cleanly observe that there Is 
Indeed a correlation between the lowest value of the SWF and the 
fracture location on the specimen. Figures 10 through 13 show similar 
correlations for various other stacking sequences for E-glass epoxy 
laminates. One should note at this point that these correlations are 
valid only for a given set of parameters. For example, if the threshold 
were different, a different fracture location might have been 
predicted. An alternate method to quantify the SWF will be discussed in 
section 2 . 2.1 where the need for selecting a threshold level and several 
other instrumentation parameters is eliminated. 

2.1.3 Correlation of SWF with Other NUT Techniques 

In the previous section, SWF demonstrated that it has the 
potential to predict failure locations on glass epoxy laminates 
nondestructi vely. To explore further the significance of the SWF 
values, various other NOT techniques such as moire interf erometry, "-ray 
radiography, ultrasonic C-scan, edge replication, etc. were applied to 
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the specimens that have been evaluated by the SWF. Of these techniques, 

A 

moire Interferometry and ultrasonic C-scan appeared to provide 
additional Insight concerning the SWF measurements. 

A total of three [O.OOjJj E-glass epoxy laminates were tested In 
this program. On each of the specimens the following tests were 
performed: 

• evaluation of the SWF along the length of the specimen, 

• ultrasonic C-scanning of the specimen, 

• calculation of the local stiffness, from the In-plane 
displacement (u) field as obtained by full field moire 
interferometry of the specimen. 

After these tests were made on the specimens, the results were analyzed 
and correlations, if any, were observed. The specimens were then loaded 
in quasi-static tension until final fracture. The fracture regions were 
checked with the locations predicted by the SWF. In this study, results 
from one specimen are presented. 

This specimen was 1.0" x 7.0" long. SWF was measured along the 
length of the specimen. Figure 7 shows a plot of the normalized SWF 
versus the position on the specimen. From Fig. 7 one can observe that 
the top end, center, and the bottom end of the specimen had low SWF 
values. Figure 14 shows the ultrasonic C-scan of the specimen. Note 
the indication of flaws at the top end, center, and the bottom end of 
the specimen. The specimen was then studied in a moire 1 nterf erometry 
set-up so that the u displacement field was obtained. Appendix A 
describes the moire interferometry technique [22]. A grating of 30,UU0 
lines/inch was glued to the specimen so that the furrows on the grating 
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were along the longitudinal axis of the specimen. The frequency of the 

virtual grating was 60,1)00 lines/inch. The effective field was about 

1.0" x 3.75". The field size Is limited by the size of the lenses. 

Since the specimen under study had a gage length of 4.0", a field of 

1.0" x 3.75" was adequate. The specimen was then fixed in the loading 
> 

frame of the moire interferometric set-up and loaded to various values 

of strain. Figure lb shows the u displacement field at 500 ye which 

corresponds to about 100 lbs. load. At this low level of strain the 

moire fringes are uniformly spaced, indicating uniform displacement 

field along the length of the specimen. Figure 16 shows the u 

displacement field at 2000 ye which corresponds to about 400 lbs. 

<* 

load. Here the moire fringes are nonuniform along the length of the 
specimen. From the displacement contours in Fig. 16 local strains were 
calculated. As the applied axial load was constant local stiffness from 
the local strains were calculated along the length of the specimen. 
Figure 17 shows a plot of the relative local stiffness and normalized 
SWF values versus the position on the specimen. From the local 
stiffness plot in Fig. 17 one can see that the top end, center, and the 
bottom end of the specimen had lower stiffnesses compared to other 
regions. Comparing Figs. 14 and 17, there appears to be a trend between 
the SWF values, ultrasonic C-scan, and the local stiffness of the 
specimen. The areas of low SWF values corresponded with the "white" 
areas of the C-scan and areas of high displacement in the moire fringe 
pattern. 

a 

The specimen, after being removed from the load frame of the moire 
interferometric set-up, was pulled in quasi-static tension until final 
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fracture. The specimen failed in the top grip. As a result, one cannot 
make meaningful correlations between fracture, the local stiffness, and 
the SWF values for this specimen. The interpretation of failure in the 
grip region is complicated because of the complex state of stress caused 
by the grips. However, it is interesting to note that the fractured 
region had low values of SWF and local stiffness (see Fig. 17). The 
other two glass epoxy laminates did fail in the vicinity of the 
predicted locations and in between the grips. Results of the specimen 
that fractured in the grip were used because the photographs of the 
moire fringes had the best contrast. Therefore, visualizing the areas 
of high and low strains is relatively easy. 

2,1.4 Variation of the SWF with Fiber Orientation 

To evaluate the effect of the variation in the SWF with fiber 
orientation, three panels: [0]g graphite epoxy, [0] 8 S-glass epoxy, and 

random fiber S-glass epoxy laminates were utilized in this study. Each 
of the panels was approximately 12 in, by 12 in. The thicknesses of the 
laminates were in the range U.U45 - 0.051 in. Figure 18 shows the 
schematic of the measurement pattern utilized. The sending transducer 
was placed in the center of the panel and the receiving transducer was 
moved around the sending transducer as a function of the fiber 
orientation. The instrumentation parameters listed in 2.1.1 were used 
here. Figure 19 shows the variation of the SWF with fiber orientation 
for a [0] 8 graphite epoxy laminate. Figure 20 shows the variation of 
the SWF with fiber orientation for a [U]y S-glass epoxy laminate. Note 
that in both cases the SWF is the maximum along the fiber direction and 
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minimum across the fibers. Figure 21 shows the variation of SWF In an 
arbitrarily chosen 0° to 90° gradient of the random fiber S-glass epoxy 
panel. Note that the SWF is approximately constant. 

2.1.5 Effect of the Number of Plies on the SWF 

Since the sending and the receiving transducers are on the same 
side of the specimen, it Is likely that the acousto-ultrasonic signal is 
representative only of the top few plies of the laminate. A simple 
experiment was performed to evaluate how many plies the SWF technique 
interrogates. Figure 22a shows a schematic of an edge replica of a 
[0,90 2 ,U] S graphite epoxy specimen. Figures 22b and 22c show the signal 
on side A and side B for this specimen. The specimen was then bent to 
introduce matrix cracks in the 90° plies towards side B as shown by the 
schematic of an edge replica in Fig, 23a. The signal corresponding to 
side A is shown in Fig. 23b. Note that there is not much change in the 
signal amplitude when compared with Fig. 22b. However, the amplitude of 
the signal corresponding to side B, Fig. 23c, decreased significantly 
when compared with that of Fig. 22c. This indicates that the SWF 
technique interrogates only the first 3-4 plies of the laminate. 

2.2 Phase 11 

In section 2.1.1, instrumentation and experimental parameters were 
listed which produced reproducible SWF values to within +10%. In 
section 2.1.2, correlations between the SWF and the fracture location on 
E-glass epoxy laminates were discussed. The correlations were found to 
hold only for the particular set of parameters listed in section 2.1.1. 
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Results are especially sensitive to the threshold level. In order to 
eliminate the need for the threshold level and most other 
instrumentation and experimental parameters, an alternative method to 
quantify the SWF was investigated. This method Is discussed in detail 
in the following section. 

2.2.1 An Alternate Method for Quantifying the SWF 

A schematic of the experimental set-up for this method is shown in 
Fig. 24. Two transducers, one acting as a transmitter, the other as a 
receiver, were mounted on the same side of the specimen as shown. The 
transducers were matched broadband transducers, 0.5" in diameter and had 
a center frequency at nominally 5 mHz. The transducers were Panametric 
model V1U9. Panametrics shear wave couplant was used between the 
transducers and the specimen. The transducers were held in a plexiglass 
fixture so that the distance between their centers was fixed at 1.0". A 
pair of elastic bands at each end of the plexiglass fixture kept the 
transducer assembly from sliding downthe specimen. A Panametrics model 
5052A was used as a pul ser/recei ver. A Panametrics model 5050AE 
preamplifier was used to amplify the received signal. A transient 
waveform recorder, Biomation model 8100, was used to digitize the 
signal. A Fast Fourier Transform algorithm was used to perform the 
spectral analysis of the signal. A typical signal and Its frequency 
spectrum for an undamaged [0 , 90g ] s graphite epoxy laminate Is shown in 
Fig. 25, Talreja [23] has noted that three classes of paremeters are 
needed to describe the distribution functions, namely location, scale 
and shape parameters. He suggested that a convenient set of parameters 
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to represent a frequency spectrum can be defined by considering the 
spectrum to be a plane figure closed on the frequency axis. The 
location parameter for the frequency spectrum can then be taken to be 
the location of the centroid of this figure. The area of this figure 
forms a suitable scale parameter and the shape parameters can be 
described In terms of the various moments of the figure about a 
convenient axis. Talreja therefore proposed that a set of shape 
parameters be defined by 

M k = 1,2,3... 

S - p- r = 2.3,4... 

r.k M r . k f c k r > k 

where 

oo 

M r = / S(f) f r d f 


where S(f) is the power spectral density, f Is the frequency, and f c Is 
the location of the centroid found by the following expression 


Several cf the moments and moment ratios can be given physical 
interpretations. For example, M 0 is the mean square value of the power 
spectral density. Al-o, 


f 


o 


■ ( 


M V 2 

V 


is the frequency of mean value crossings with positive slopes and 



is the frequency of maxima in the time domain [24], In this study the 
root mean square of the power spectral density, (M )^ p was used as the 
parameter for the alternate method to quantify SWF. In applying this 
technique the need for the threshold level and several other 
instrumentation parameters is eliminated. The only instrumentation 
parameter used was the gate width. The gate width was chosen so that 
only the initial part of the signal was received by the transient wave 
form recorder. 

2.2.2 Correlation of the Root Mean Square of the Frequency Spectrum 
with the Laminate Stiffness During Fatigue 
A schematic diagram of the experimental set-up is shown in 
Fig. 24. The specimens tested were T-300/5208 graphite epoxy laminates 
measuring 1.0 in. by 8.0 in. long. Three stacking sequences were 
Investigated: [0,90 2 ] s , [0,+45] s and [0,9U,+4t>] s . The ends of the 

specimens, approximately 1.5 in., were coated with a thin uniform layer 
of MS907 epoxy manufactured by Mi 1 ler-Stephenson Chemical Company. This 
coating provided a uniform gripping surface. TV'n thicknesses of 320 
grit emery paper were inserted between the MTS grips and the specimen 
during fatigue loading to improve the gripping efficiency. The cyclic 
loading was done on an MTS servo-hydraulic, closed-loop testing frame. 
Stiffness was measured by means of an extensometer attached to the 
center portion of one side of the specimen. The extensometer gage 
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length was nominally 2 In. The extensometer knife edges sat In V-shaped 
grooves machined Into aluminum tabs. These tabs were bonded to the 
specimen by means of RTV silicone rubber cement manufactured by 3-M 
Company. On the other side of the specimen, the transducer assemply 
discussed in section 2.2.1 was attached. The transducer assembly was 
positioned so that It was at the center of the specimen directly 
opposite the extensometer. 

The fatigue tests were conducted at a constant stress ratio, R“0.1, 
In a sinusoidal form at a cyclic frequency of 10 Hz In a tension-tension 
mode. The maximum stress amplitudes used for the testing [2] were as 
follows: 

C0.90 2 a s . o nax • 0.70 S u|t 

ro.*45] s . 0max • 0.72 S uIt 

C0.9O.t4S],, " 0.62 S u , t 

where Is the ultimate tensile strength. The pulser transmits 
ultrasonic pulses through the sending transducer. The ultrasonic signal 
was picked up by the receiving trandsucer. The signal was sent through 
i preamplifier and then an amplifier. A filtering circuit was used to 
filter out frequencies less than 100 kHz. This removed any noise 
associated with the MTS hydraulic system. An oscl 1 loscope was used to 
display the signal. The signal was then fed Into a counter where the 
counts above a certain threshold were calculated (conventional SWF). 

The signal was also captured and digitized using a transient waveform 
recorder. The signal was stored on a floppy disk of a Z-80 based 
microprocessor. A Fast Fourier Transform algorithm was used to generate 
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the frequency spectrum of the signal. Root mean square of the frequency 
spectrum {also called power spectral density) was calculated from the 
spectrum. 

Results from a typical [O.yOgjg graphite epoxy specimen are 
presented next. The letter N shall be used to denote the number of 
cycles, and P, the load. The specimen was Inserted Into the grips of 
the MTS machine and all the necessary attachments were made to obtain 
the signal and the stiffness,. Figur® 26 shows the signal and the 
frequency spectrum at N=U and P°0. Note the presence of two harmonic 
frequencies, one centered approximately at 300 kHz and the other 
centered approximately at 7b0 kHz. Figure 27 shows the signal and the 
frequency spectrum for N=V2» P = mean load, 1248 lbs. Note the reduction 
in the higher harmonic frequencies. Figures 28 through 34 show the 
signals and their corresponding frequency spectra for various numbers of 
cycles and at P-mean load, throughout the fatigue life of the 
specimen. In all cases signals and stiffnesses were recorded at mean 
load with the cycling stopped. Root mean square values of the spectra 
were then calculated. Table 1 shows the root mean square values and the 
stiffnesses for this specimen. Figure 3b shows the stiffness and the 
root mean square value versus the number of cycles. Note the 
correlation between the stiffness and the root mean square values of the 
frequency spectra. Figure 36 shows a plot of fp versus strain for 
[G,902 ] s laminate. Figure 37 shows a plot of fc and fo versus strain 
for [0 , 902 H s laminate. Figure 38 shows penetrant-enhanced X-ray 
radiographs of [O.l^jg graphite epoxy laminate at 60,000 cycles and 
100,000 cycles. 
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Figures 39 through 4b show the signals and the corresponding 
frequency spectra for various number of cycles for a [U,t4b] s graphite 
epoxy laminate. Figure 46 shows the plot of stiffness and the root mean 
square values versus the number of cycles for [0,i4bL laminate. Note 
the correlation between the stiffness and the root mean square values of 
the frequency spectra. 

Figures 47 through 51 show the signals and the frequency spectra 
for various number of cycles for a [0,90,±45] s laminate. Figure 52 
shows the plot of stiffness and the root mean square values versus the 
number of cycles for [0,9(J,t45] s laminate. Note the correlation between 
the stiffness and the root mean square values of the frequency spectra. 

2.3 Phase III 

In section 2.2.2, correlations between the root mean square of the 
frequency spectrum and the stiffness reduction during fatigue were 
shown. The root mean square values of each of the spectra were 
calculated manually. The number of signals and their spectra would get 
large for long-term fatigue tests. To decrease the computation time 
necessary to calculate the root mean square values, a computer code was 
written. The quantities that the program calculates are 
(M 0 ) /2 , (M^/M^ 2 . Section 2.2.1 describes these 

ratios in detai 1 . 
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2.3.1 Computer Code to Calculate Various Moments of the Frequency 
Spectrum 

A pu-:er sends ultrasonic pulses through tne specimen by means of a 
sending transducer. The signal received by the receiving transducer Is 
preampl If led, amplified, and after appropriate signal conditioning, Is 
sent to a transient waveform recorder where the signal Is digitized Into 
2048 points and stored. The data Is then sent through a parallel I/O 
that Is 8 bits wide, to a microprocessor. From the microprocessor the 
data Is sent through a serial I/O and a modem to an IBM 37U mainframe 
computer for storage on disk. The entire process of sending a pulse 
through the specimen to storing the signal on disk takes about 50 
seconds. An IMSL [25] subroutine wcs used to calculate the Fast Fourier 
Transform of the signal. To calculate the various moments of the 
frequency spectrum, a computer code was developed. Appendix B has the 
computer program. 

2.3.2 Flow Chart of the Data Acquisition, FFT, and the Moments 

Evaluation System. 

The flow chart for the system to capture the signal, transformation 
to the frequency domain, and calculation of the moments of the frequency 
spectrum is shown in Fig. 53. 


3. DISCUSSION AND SUMMARY 

This study has Investigated the use of stress wave factor as a NDT 
technique for thin composite laminates. For the SWF to be a useful NUT 
technique, it must possess the following requirements: 

• provide reproducible results 

• should have very few experimental and instrumentation 
parameters 

« should be easy to operate 

• should be operator insensitive 

• the results achieved by the technique should ha'/e the 
potential to be quantified 

• the results should provide insight Into the nature of damage 
development 

• the technique should be able to yield a parameter that Is 
indicative of the total damage between the transducers 

• this parameter should be able to correlate/predict mechanical 
behavior 

• the technique should have the capability of automatic 
Inspection 

From the results presented In the previous chapter, most of the 
above requirements appear to have been mat. The technique provides 
reproducible results, is easy to operate, is operator Insensitive, 
provides quantitative parameters, and provides some Insight into the 
nature of damage development in composite laminates. Additional work 
needs to be done before one can be sble to relate the quantitative 
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parameters obtained bj this technique with the mechanical response of 
the material. However, a description of the nature of wave propagation 
Is difficult because of the anisotropy of the composite laminates and 
the complex damage state that develops as a result of mechanical 
loading. When a mechanical wave propagates through a material it 
undergoes reflections, mode conversions, several attenuation mechanisms, 
dispersion, etc. Any variation in the material condition will generally 
cause a change in one or more of the above areas. A precise description 
of the material 's condition is not possible using these concepts. If 
one were to have a precise description of the material 's condition then 
the description will at best facilitate determination of the stiffness 
since the strength or lif* depend on the nature and history of the 
applied loads. In a large structure, determining the material condition 
everywhere is impractical. It is therefore desirable to determine where 
the material condition Is "poor". Knowing this, attention can then be 
directed to determine the exact condition. It appears that the SWF 
technique as applied in this study possesses a tremendous potential for 
doing just this type of inspection. 

The results presented in the previous chapter will be discussed in 
three sections. Section 3.1 deals with the initial phases of the 
investigation which includes reproducibi lity, failure location 
prediction, correlation with other NDT techniques, effect of fiber 
orientation on the SWF, effect of the laminate thickness on the SWF, and 
the nature of wave propagation. Section 3.2 discusses the frequency 
analysis and the correlation of the root mean square values to the 
laminate stiffness and damage development during cyclic loading. 
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Section 3.3 discusses the usefulness of the SWF as an NDT technique, 
advantages and limitations of the SWF. 

3.1 Initial Investigation 

For any NDT technique to be useful, It must be able to provide 
reproducible results. Using the AET model 206AU, several carefully 
monitored experiments were conducted on E-glass epoxy laminates. The 
outcome of these experiments was a set of parameters that yielded 
reproducible SWF values to within ±10% (see section 2.1.1), The 
variation of the SWF values by as much as 10% is due to the inability to 
place the transducers in the same location on the specimen, variation in 
the amount of couplant used, effects of temperature and moisture, and 
misalignment of the transducers relative to the specimen. 

When virgin E-glass epoxy laminates wre evaluated by the SWF and 
then pulled In quasi-static tension, the specimens failed in the regions 
where the SWF was the lowest. This behavior suggests that areas of low 
SWF would be indicative of areas with poor stress wave energy 
transmission, thereby dynamic strain energy is likely to concentrate in 
such areas and promote failure. Local stiffnesses calculated from full 
field moire Interferometry correlated well with SWF for E-glass 
laminates. This suggests that the local 1 nhomogeneitl es present in the 
laminate that cause variations in mechanical properties also Influence 
the stress wave propagation and hence the SWF. Ultrasonic C-scans of E- 
glass laminates were found to complement the SWF results. It Is 
interesting to note that in the limiting case where a single transducer 
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Is used as a sender and a receiver the SWF technique becomes the 
conventional pulse-echo A-scan. 

The SWF was found to be sensitive to the fiber direction of the 
composite laminate. When unidirectional, 8-ply graphite epoxy and 
S-glass epoxy laminates were evaluated by the SWF, the SWF values were 
found to be the largest along the fiber direction and the smallest 
perpendicular to the fiber direction. This behavior suggests that the 
0° plies transmit the stress wave energy more efficiently compared with 
90° plies. It appears, therefore, that fiber orientation Is critical 
not only to sustaining static loading but also to dynamic stress wave 
transmission. 

A sequence of tests performed on [0,90?,Q] S graphite epoxy 
laminates have shown that the SWF obtains information mostly from the 
first four plies of the laminate. 

3.2 Frequency Analysis 

To eliminate the need for several experimental and instrumentation 
parameters, livted in section 2,1.2, the root mean square of the 
frequency spectrum of the received signal was investigated as an 
alternative method for quantifying the SWF, rather than counts above an 
arbitrary threshold level. This method has better sensitivity and is 
more practical to use than the conventional SWF. Figure 31 shows a 
signal and its spectrum for a CO, 90^2 graphite epoxy laminate that has 

„N 

been cycled to 30,000 cycles. The SWF Is 787 and the (-£—) Is 0.68. 
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Figure 32 shows the signal and its spectrum for the same specimen at 
95,000 cycles. Note that there is little change in the signal. The SWF 

H N '/2 

is 790 while the (-— ) is 0.61. It has been shown [1] that for this 

M o 

type of laminate and test conditions there is a significant growth of 
micro-delamination from around 30.0UU cycles to 10U,0UU cycles. The 
root mean square of the frequency spectrum is therefore more sensitive 
to damage than the conventional SWF. It is also more useful because a 
great deal of Information is obtained about the signal from the 
frequency spectrum. One can obtain the root mean square value, the 
approximate frequency of mean zero value crossings with positive slopes, 
the approximate frequency of the maxima In the time domain and the 
centroid of the frequency spectrum (see section 2.2.1). The results 
presented in section 2.2.2 show that the root mean square values of the 
frequency spectrum correlate extremely well with the reduction in 
laminate stiffness during cyclic loading. Three different stacking 
sequences of graphite epoxy laminates were tested. Results from each of 
these stacking sequences is discussed next. 

[0,90 2 ] s Laminate 

Figure 26 shows the signal and its frequency spectrum of the 
specimen at N=U and P S U. Note the presence of two harmonic frequencies, 
one centered at approximately 25U kHz and the other centered at 
approximately 7b(J kHz. Figure 21 corresponds to N=V 2 and P^mean 
load s 1238 lbs. The amplitudes of the second harmonic frequencies are 



28 


reduced. It has been shown [2] for the test conditions used, transverse 
cracks in the yd' plies appear at this load level. Therefore the 
transverse cracks appear to cause a reduction in the high frequencies 
associated with the wave. Figure 28 corresponds to N=2000 cycles. The 
frequency spectrum of the second harmonic frequencies has disappeared. 

It has been shown [2] that the characteristic damage state is reached 
around a few thousand cycles. Therefore the presence of transverse 
cracking has resulted in a complete reduction in high frequencies 
associated with the wave. Note also that the frequency spectrum of 
first harmonic frequencies has become skewed and its amplitude has 
decreased. The skewness of the spectrum increases until about 35, '"JO 
cycles. Jamison [2] has shown that around this level of cycling the 
micro-delaminations that have been growing reach a fairly widespread 
distribution in the entire laminate. Around N=35,OUO cycles the 
frequency spectrum becomes Gaussian again. As cycling continues, micro- 
delaminations enlarge gradually and the amplitude of the spectrum 
decreases gradually, Figs. 29 to 31, At N=95,000 cycles, Fig. 32, there 
is a significant drop in the amplitude of the frequency spectrum. It 
has been shown [2] that around 100,000 cycles, micro-delaminations tend 
to join together to form internal macro-del ami nations. Figure 38 shows 
penetrant-enhanced X-ray radiographs of a 1^0 ,90g] s graphite epoxy 
laminate used in this study. Figure 38a corresponds to N=70,000 
cycles. Note small fuzzy areas around where the cracks intersect. 

These are micro-delaminations. Figure 38b corresponds to N=100,000 
cycles, after a drop has been noticed both in the stiffness and the root 
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mean square values. Note the small .Teas of micro-delaminations in the 
interior of the specimen joined together to form macro-delaminations. 

At N=140,000 cycles. Fig. 33, the amplitude of the spectra 
decreases further. After this the amplitude remains nearly constant, 
Fig. 34, with cycling reflecting the Stage II behavior of this 
particular laminate [2], The root mean square values of these frequency 
spectra correlated well with the reduction in stiffness load cycles 
because of damage, Fig. 35. Note that the reduction in the root mean 
square values is about 6U% whereas the reduction in stiffness because of 
damage is only about 13%. Therefore the root mean square of the 
frequency spectra exhibits a larger decrease than the stiffness. 

For [0,90 2 ] s graphite epoxy laminates transverse cracks occur at 
0.4% strain. When f Q , f c , and f p (see section 2.2.1) were plotted 
versus strain, there was a sudden decrease in these values at 0.4% 
strain, Figs. 36 and 37. This suggests that the parameters f Q , f c , and 
fp, which are ratios of various moments of the frequency spectrum could 
also be used to monitor damage development in composite laminates. 

[0,±4b] s Laminate 

Figure 39 shows the signal and the frequency spectrum for a 
[U,±45] s graphite epoxy laminate at N=0 and P=0. The frequency spectrum 
for this specimen is not as well defined as for the [G,90 2 ] s . The 
frequencies are spread between 100 kHz and 1 MHz. Figure 40 corresponds 
to N= 1/2 and P=mean load=1345 lbs. There is no significant change in 
the spectrum. The strain levels are not sufficient at this load to 
cause matrix cracking in the ±45° plies. Figure 41 corresponds to 
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N=1000 cycles. There Is a small decrease in the amplitude of the 
spectrum over the entire frequency range. Figures 42 through 45 
correspond to N=21,00U cycles, 35,000 cycles, 45,000 cycles and 95,000 
cycles. The amplitude of the spectrum gradually decreased with 
increasing cycles. It has been reported [2] that the damage modes in 
these laminates are restricted to matrix cracking in the ±45° plies and 
longitudinal splitting in the 0° plies. There is no significant 
delamination. This is believed to be one of the reasons why a sudden 
drop in the (M 0 )^ a nd stiffness values with cycling as in the case of 
[0,90 2 ] s laminates does not occur. Figure 4b shows that the root mean 
square values of the spectra correlated quite well with the stiffness 
reduction with cycles. Note that the exhibits a larger decrease 

than stiffness. 

[Q,90,±45] s Laminate 

Figure 47 shows the signal and its frequency spectrum at N=0 and 
P=0. The spectrum appears to be somewhat between that of a [0,9U2] S and 
a [0,±45] $ spectra. Figure 48 corresponds to N=V 2 and P=mean load=l254 
lbs. There Is no significant change in the spectrum. Figure 49 
corresponds to N=2U,0U0 cycles. There is a small decrease in the 
amplitude of the frequency spectrum. At N=30,000 cycles, Fig. 50, the 
amplitude drops by a considerable amount. It has been reported [2] that 
around 25,000 cycles, delaminations grow slowly from the free edges into 
the specimen. Visual confirmation of edge delamination was made at 
N=50,000 cycles. Fig. 51, for the specimen being tested. Figure 52 
shows that the root mean square values of the spectra correlate quite 
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well with reduction In stiffness with load cycles, because of damage. 
Note that the (M q )^ 2 exhibits a larger decrease than stiffness. 

Therefore, in all three laminate types tested, there was an 
excellent correlation between the root mean square values, stiffness, 
and damage development in controlled cyclic loading. In all cases, 
the (M 0 ) 1/2 exhibited larger decreases than the stiffness. A computer 
program to calculate (Mq)^ 2 , fp, f c , and f 0 was developed to save 
computation time. The signal received by the receiving transducer, 
after appropriate signal conditioning, is sent through a transient wave 
form recorder and then to the computer. The computer program calculates 
(M 0 ) /2 , fp, f c , and for the spectrum of each signal (see Fig. b3). 

3.3 Usefulness of the SWF as an NOT Technique 

From the results and discussion reported thus far. It is clear that 
there Is Indeed a correlation between ( M Q } , stiffness, and damage 
in thin composite laminates. The SWF technique is particularly useful 
because it yields a single parameter, such as ( M q ) , to characterize 
the damage, and the effects of various forms of damage between the 
transducers in Integrated, very much like the stiffness. From an 
engineering standpoint this feature has tremendous significance. The 
parameter so obtained could be used in a damage parameter-stiffness-life 
model [26] to predict the mechanical response of the material because of 
the presence of damage. 
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Advantages of the SWF Technique 

• This technique yields quantitative parameters, 

like (M , which can be correlated with the stiffness of 
composite materials. 

• (rg 1 ^ exhibited larger decreases than stiffness for cyclically 
loaded composites. 

• It requires access to only one side of the structure, which is 
of great practical Importance. 

• It has the potential to be used as an integral NDT system with 
a structure to monitor damage development. 

Limitations of the SWF Technique 

• This technique cannot provide a precise geometric description 
of damage In the material, 

• The strength of the signal changes with the distance between 
transducers. 

• A unique relationship between strength and the parameter 
obtained by SWF is not possible. 



33 


4. CONCLUSIONS 

The conclusions arrived at from this study are listed below. 

• The SWF has excel lents potential as an NDT technique for thin 
composite laminates. 

• The results of other NOT techniques such as moire 
Interferometry and ultrasonic Oscan support the SWF results. 

• SWF Is the highest along the fiber direction and the lowest 
across the fibers. 

• SWF appeared to interrogate the top 3 to 4 plies of the 
laminate for the test conditions used in this study. 

• An alternative "SWF" method, which makes use of the root mean 
square values of the frequency spectrum of the received signal 
was developed. 

• This alternate method does not need a threshold level and 
hence Is more practical, 

• The root mean square of the frequency spectrum correlated with 

the reduction In the laminate stiffness during cyclic loading 
for the three stacking sequences of the laminates were 
tested: [0,90 2 ] s , [0,±45] s , and [0,90,±45] s . 

• The root mean square values exhibited larger decreases than 
stiffness, during cyclic loading of composite laminates. 

• The root mean square values have potential to provide Insight 
Into the nature of damage development. 
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Recommendations for Future Work 

• Perform basic experiments to understand the nature of wave 
propagation In composite laminates. 

• Perform systematic experiments to correlate the exact nature 
of damage (as seen by X-ray radiography, edge replication, 
etc.) with the various moments of the frequency spectrum. 
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Table | 


N 

E 

Specimen [0,90 2 ] 

E/E 0 

s Gr/Ep 

M 

N/N 0 

0 

7.34xl0 3 (»E o ) 

1.00 

493.25xl0” 5 (=M o ) 

1.00 

l ll 

7.12xl0 3 

0.97 

428.75xl0" 5 

0.93 

2,000 

6.97xl0 3 

0.95 

307.84xl0“ 5 

0.79 

6,000 

6.90xl0 3 

0.94 

262.85x10“® 

0.73 

16,000 

6.83xl0 3 

0.93 

248.65x10“® 

0.71 

30,000 

6.83xl0 3 

0.93 

202.03x10“® 

0.64 

95,000 

6.53xl0 3 

0.89 

195.77x10“® 

0.63 

140,000 

6.46xl0 3 

0.88 

82.91x10“® 

0.41 

270,000 

6.38xl0 3 

0.87 

67,52x10“® 

0.37 


All (Measurements made with cycling stopped and P=mean load«1248 lbs. 
Units for E are ksl; units for M are volts/sec. 
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Fig. 1 


Schematic of the SWF set-up 
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Fig. 5. Narrow gate width at two different threshold levels yielding 
approximately the same number of counts for two substantially 
different signals. 
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g. 6. Schematic showing the location of the SWF measurements made on 
the specimen. 
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Plot of the normalized SWF vs. position on spec 
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Fig. 10 


Correlation of SWF with specimen failure location 
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Fig. 11. Correlation of SWF with specimen failure location. 
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Fig. 14. Ultrasonic C-scan of a virgin [0,90^] E-glass epoxy 

specimen. Note Indication of flaws at the top end, center, 
and the bottom end of the specimen. 
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* 

Fig. 15. Moire fringe pattern (u displacement field) of a [0 , 90 3 ] s 
E-glass epoxy specimen at 500 Note the uniform 

displacement field along the length of the specimen. 
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Fig. 15. Moire fringe pattern (u displacement field) at a [0 , 9U 3 3 S 
E-glass epoxy specimen at 2000 ,ic. Note that the top end, 
center, and the bottom end of the specimens have higher 
displacements compared to other regions. 
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Plot comparing normalized SWF and local stiffnesses vs 
position on specimen. 
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Fig. 13. Schematic of the set-up to evaluate the effect of 
orientation on the SWF, e = sender, r = receiver 
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Fig. 20. Variation of the SWF with fiber orientation. 
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Fig. 21. Variation of the SWF in an arbitrarily chosen 0° to 90° 
quadrant of a random fiber S-glass epoxy panel. 
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Fig. 22b. ^^fto-ultrasonic signal on side A of an undamaged [0,90 2 ,0^s 
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Fig. 22c. Acousto-.Ur.sonU signal on side 8 of an undamaged t0.9O £ .03« 
specimen . 
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Fig. 23b. Acousto-ultrasonic signal on side A of a damaged [0,90 2 ,0) s 
specimen. 
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Fig. 23c. Acousto-ultrasonic signal on side B of a damaged [0,90 2 ,0] s 
specimen. 
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Fig. 24. Schematic of the experimental set-up to calculate the Fourier 
Transform of the signal. 
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Fig, 26. The signal and Its spectrum at N=0, °*0 for a [0,902] s 
graphite epoxy specimen. 
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Fig. 28. The signal and its spectrum at N=2,000, P=1248 lbs for a 
[0,90 2 ] s graphite epoxy specimen. 
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Fig. 29. 


The signal and its spectrum at N=6,QUU, P=1248 lbs for a 
[0,902 ] s graphite epoxy specimen. 
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Fig. 30. The signal and Its spectrum at N=16,000, P=1248 lbs for a 
[O.Wgig graphite epoxy specimen. 
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Fig. 31. The signal and its spectrum at N=30,000, P= 1243 lbs for a 
[0 , 902 3 S graphite epoxy specimen. 
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Fig. 33. The signal and its spectrum at N®140,000, P=1248 lbs for a 
[0,9t)2] g graphite epoxy specimen. 
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Laminate stiffness (solid line) and root mean square of the 
frequency spectrum (dotted line) for fatigue test of a 
[0,90 2 ] s graphite epoxy specimen. 
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Fig. 36. Frequency of maxima In time doman, fp vs. strain for a static 
test of a [0,90 £!) s graphite epoxy specimen. 
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Fig. 37. Frequency of the mean value crossings, fo, and the location of 
the centroid fc, vs. strain for a static test of a [0,90 2 ] s 
graphite epoxy specimen. 
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Fig. 38. Penetrant-enhanced radiograph showing (a) mlcro-delamlnations 
and (b) macro-delaminations In a [0,90?] s graphite epoxy 
specimen that has been cyclically loaded. 
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Fig. 40. 


The signal and its spectrum at N=^j P= 1 345 lbs for a [0,±4S] s 
graphite epoxy specimen. 
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The signal and its spectrum at N=1000, 1 345 lbs for a 

[0,±4S] s graphite epoxy specimen. 
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Fig. 45. The signal and its spectrum at N=95,000, P=1 345 lbs for a 
[0,±45] s graphite epoxy specimen. 




tC/M 


85 


oRiGtf'iAL 
OF 'POOR Ql'.V .V v 



N, kc 


Fig. 46. Laminate stiffness (solid line) and root mean square of the 
frequency spectrum (dotted line) for fatigue test of a 
[0 , ±45 3 S graphite epoxy specimen. 
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Fig, 47, The signal and its spectrum at N=0, P=0 for a [Q,90,i45] s 
graphite epoxy specimen. 
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Fig. 48. The signal and Its spectrum at N =*& P=1254 lbs for a 
[0,90,±45] s graphite epoxy specimen. 
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Fig. 49. The signal and its spectrum at N=20,000, P=1254 lbs fir a 
[0,90,±45^) s graphite epoxy specimen. 
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Fig. 52. Laminate stiffness (solid line) and root mean square of the 
frequency spectrum (dotted line) for fatigue test of a 
[0,90,±45] s graphite epoxy specimen. 
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APPENDIX A 

Molrg interferometry depends upon diffraction of light as wiiH as 
interference [ ]. Moire fringes are obtained by using diffraction 
gratings. A grating is a surface with regularly spaced bars or furrows. 
The distance between two consecutive bars is called pitch, g. Frequency 
f of a grating is the number of bars per unit length 



A grating divides every incident wavetrain into a multiplicity of 
wavetrains of smaller intensities; and it causes these wavetrains to 
emerge in certain preferred directions. A parallel beam incident at a 
particular ongle on the grating divides it into a series of beams that 
emerge at preferred angles. These beams are called diffraction orders 
and are numbered in sequence beginning with the zero order, which is the 
mirror reflection of the incident beam. The angle between the neigh- 
boring < ./fraction orders is small for a coarse grating and it is large 


laser 



Figure a. - Schematic of moire Interferometry setun. 
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for a fine grating. 

The grating equation defines the angles of diffraction, viz.. 

Sin 0 * mX + Sin a 

where m defines the diffraction order of the beam under consideration. 

Coarse moire Is produced by the superposition of two coarse gratings 
with frequencies In the range of about 1 to 40 llnes/ran. One of them is 
an active grating (glued to the surface of the specimen) and the other a 
reference grating (virtual). The virtual grating is formed by the 
intersection of the wavetrains from the reflective active grating and a 
mirror beside the a:t1ve grating. Fig. a shows the moird interferometry 
set up. 

Molrd fringes are the locu*- of points of cnnstant displacements, 
specifically the In-plane displacement component in the direction 
perpendicular to the lines of the reference grating [ 22 \, Molrd fringe 
order, N, denotes the number of cycles of intensity fluctuation experi- 
enced at any point as the displacement changes from zero to its final 
value. 

u * g N ; x = longitudinal direction 

'X 

v ~ g Nyj y = transverse direction 

--u , v are displacements in the longitudinal and transverse directions 
--g a pitch of the reference grating 

— N x , N v are the fringe orders when lines on the reference grating are 
perpendicular to the longitudinal and the transverse directions, re- 
spectively. 

From the displacements we can get strains 
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. . du 
e x 37 

_ „ dv 

y 37 

r . + dv 

yy dy dx 

where t ard e are the normal strains and y is the shear strain, 
x y (y 

Stresses can be determined using stress-strain relationships. Since 
the stains and the associated displacements are very small, gratings of 
fine pitch and high frequency are required C 22 l • 

A grating is made by exposing a high-resolution photographic plate 
to two intersecting beams of coherent light. Frequency of the grating 
is controlled by the angle of intersection according to 


Sin a m — f 

Ct 

where a = incident angle 

X a wavelength of light used 
f = frequency of grating 

The two intersecting beams give rise to a diffraction pattern of light 
and dark bands. When the plate is developed silver grains remain in the 
exposed zones, while the silver is leached out in the unexposed zones. 
The gelatin matrix shrinks upon drying, but since it is partially 
restrained by the silver, shrinkage is greatest in the unexrused zones. 
The result is a plate with a furrowed surface which can be used as a 
grating. The grating is mirrorized with aluminum to improve its reflec- 
tive qualities. The grating is then transferred and attached to ’he 
specimen with an adhesive. 
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APPENDIX B 


Computer Program 

//B122SPUT JOB 11992, PAS I, REC10N=250K, TIME=( 0.59) , M5GLEV2L= (0,0) 

/•PRIORITY IDLE 

/* JOBPARM LINES=100, CARDS-100 

//STEP 1 EXEC FORTXCGV 

//FORT. SYSPRINT DD DUMMY 

//FORT. SY1 IN DD » 

C$ JOB ANIL,PACE=20 

IMPLICIT REAL*8(A-H,0-Z) 

INTECER N, I WK (11) , I A(2048 ) 

REAL AI(2048),WK(1), AMP ( 1025 ) 

COMPLEX X( 1025 ) 

C 

C V0L= INPUT VOLTAGE RANGE 

C 

DATA PI , V0L,DELX/3 . 141592654D0, 2 , 0D0, 0. 05D0/ 

N=2048 

NNH=N/2+l 

TOTAL=N*DELX 

READ 100,(IA(I),I=1,N) 

DO 10 KK=1,N 
AI (KK)=0, 0 
AI(KK)=IA(KK> 

IF(AI(KK).CT.127.0)AI(KK)=AI(KK)-256.0 
10 A I (KK)=VOL # AI (KK)/128.0 
CIO A I (KK)=VOL*(AI(KK, -128.0)/128.0 
PRINT 202 
DX=O.ODO 
DO 20 IJ=1,N 
PRINT 110,DX,AI(IJ) 

20 DX=DX+DELX 
C 

CALL FFTRC(AI,N,X, IWK.WK) 

C 

C PRINT 110, (X(I ), 1=1, NNH) 

C 

PRINT 201 

DO 30 1=1, NNH 

AMP( I ) =0 . 0 

FN=2 . 0*PI * I/TOTAL 

AMP ( I ) =CABS ( X ( I ) ) 

30 PRINT 11C, FN, AMF{ I ) 

C THE FOLLOWING IS FOR CALCULATING THE MOMENTS 

NM=NNH-1 

FNT=2.0*PI*NNH/TOTAL 

AM0=0 . 0 

AM1=0, 0 

AK2=0.0 

AM4=0 . 0 

DO 40 1=1 ,NM 

AMO=AMO+ ( AMP ( I +1 ) -AMPII ) ) * ENT/2 . 0 

AM1=AM1+ ( AMP (1*1) -AMP { I ) ) *FNT/2 .0 * ((I-1)*FNT + FNT/2,0) 
AM2=AM2+(AMP< I«-1)-AMP(I) )*FNT/2.0 * ((I-1)*FNT ♦ FNT/2.0)**2 
40 AM4=AM4+{ AMP( I ♦ 1 ) -AMP( I ) ) *FNT/2 . 0 « ((I-1)*FNT * FNT/2.0)**4 

A=DSQRT( DABS (AMO)) 

B=AM1/AM0 

C=DSQRT(DA!iS(AM2/AM0) , 

D=DSQRT ( DABS ( AM4/AM2 ) ) 

PRINT 130, A, B , C, D 
C 

100 FORMAT(I3, 1714) 

110 FORMAT (2E 15. 5) 

120 F0RMAT( E1S . 5 ) 

13J F0RMAT(/, ’ SQUARE ROOT OF M0= ’ , E1S . S , /, 1 Ml OVER M0= ’ , E15 . S , /, 

* 1 SQUARE ROOT OF M2 OVER M0= ’ , E15 . 5 , /, 

♦ ’ SQUARE ROOT OF M4 OVER M2= ' , E15 . 5 , /) 

201 FORMAT(/, 'FREQUENCY VS AMPLITUDE' ,/) 

202 FORMAT{/, 'TIME VS VOLTAGE' ,/) 

STOP 
END 

//GO. SYS IN 
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